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Abstract 

An objective account of the action of a Stern-Gerlach apparatus on spin-1 
particles is given, using the Dirac equation. This generalizes earlier work on a 
causal interpretation of the Pauli equation to the relativistic domain, leading 
to a more natural choice for the current in the model. 


1 Introduction 

Dewdney et al. [1, 2] have described the action of an inhomogenous magnetic field on 
an uncharged spin-| particle, as in a Stern-Gerlach type spin measurement. Their 
aim was to demonstrate that there is a causal interpretation of quant um mech¬ 
anics which is not only consistent with the equations of quantum mechanics, but 
which quantitatively reproduces the results of measurements made on ensembles of 
identically prepared systems. According to such an interpretation we may assign 
well-defined trajectories and spin vectors to particles and study how these interact 
with the apparatus during the course of a measurement. For a spin measurement, 
the conventional collapse of the wavefunction onto eigenstates of the spin operators 
along the field direction is replaced by a continuous evolution of the particle’s spin 
vector towards either alignment or anti-alignment with the field direction. Dewd¬ 
ney et al. interpret this alignment as the consequence of ‘quantum torques’ which 
act in addition to the usual spin precession in the magnetic field. These quantum 
torques depend on the wavefunction and hence are sensitive to the whole experi¬ 
mental context. 

In this letter we provide a fully relativistic account of a similar spin measurement. 
It is well known that the current employed by Dewdney et al. is inconsistent with 
that obtained from Dirac theory in the non-relativistic limit, the two differing by a 
term in the curl of the spin vector [3, 4]. However, the qualitative results of Dewd¬ 
ney et al. for the behaviour of the streamlines and the spin vector coincide with our 
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predictions in the non-relativistic case. Our results demonstrate that the determin¬ 
istic evolution of the wavefunction suffices to account for the discrete outcomes of a 
quantum spin measurement. We do not need the relativistic analogue of Bohm/de 
Broglie theory in order to accept the validity of these results. Besides dealing with a 
well-defined current, a relativistic treatment is necessary as a foundation for future 
work extending these ideas to correlated spin measurements made over spacelike 
intervals, as in EPR-type experiments, so as to take issues of locality properly into 
account. The treatment presented here is, at the first quantised level, appropriate 
to scattering from arbitrarily large magnetic field gradients. However we shall see 
that pair production becomes significant in such cases, demanding ultimately a field- 
theoretic treatment. With an appropriate change in the electromagnetic interaction 
term, the present treatment could be extended to model the motion of an electron 
in an inhomogeneous external magnetic field. This would shed light on Bohr’s state¬ 
ment that no measurement of the magnetic moment of the electron could be made, 
by simply observing the beam splitting on passage through such a field (see e.g. 
appendix to [5]). 

In the extension of Bohm’s original ideas to include the concept of spin [6], the 
Pauli spinor is parameterised as 


|$) = Re ^/ 2 
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where {0, <f>, ip} are Euler angles describing the state of rotation of the particle. This 
spinor is assumed to evolve according to the Pauli equation, 


iW, |$) = -|^v 2 |$) 


( 1 . 2 ) 


where p is the anomalous magnetic moment of the particle. The particle velocity v 
is defined by J = pv, where p = ($|$) and the current J has components J k given 
by 



A = w#>) - 


(1.3) 


Note that no spatial integration is implied by the use of bra-ket notation. The 
spin-vector s has components Sk given by 
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where the {d^} are the Pauli spin operators. The Pauli equation can now be used 
lo determine the equation of motion for the spin vector along a chosen streamline. 
/The resulting equation differs from the classical equation by the inclusion of the 
/ ‘quantum torque 5 which is a non-local term involving spatial derivatives of the spin 
current ps. 
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Using these definitions of s and v, Dewdney et ah calculate the integral curves 
of the velocity vector (streamlines) and the Euler angles {0, <^>, ip} as functions of 
spatial position and time for a Gaussian wavepacket which is localised along the 
field direction. They find that the wavefunction evolves in just the way required to 
give either alignment or anti-alignment with the field direction. 

Our treatment of the relativistic problem uses the Dirac equation expressed in 
the spacetime algebra (STA) formalism of Hestenes [7], which we briefly review in 
the next section before applying it to the spin measurement in Sections 3 and 4 . 
We believe that the STA is not only a powerful computational tool for discussing 
quantum mechanics but one which brings genuine new insight into the geometric 
structure of Dirac theory, implying far-reaching consequences for the quantum the¬ 
ory of measurement. Of course, the results presented here are independent of the 
representation in which the underlying theory is expressed, but we believe their in¬ 
terpretation is most transparent when considered in the STA. We use natural units 
(c = 1 , h = 1 ) except when presenting the results of our simulations. 

2 Dirac Theory in the STA 

The spacetime algebra [7, 8 ] is a real, geometric (Clifford) algebra [9] on Minkowski 
spacetime. Vectors a, b are equipped with a single product that is associative (and 
distributive over addition), which we denote ab. It is convenient to introduce four 
orthonormal basis vectors { 7 ^}, // = 0 , 1 , 2 , 3 , satisfying 

7m'7^ = |(7m7m + 7^ 7m) = = diag (+-), (2.1) 

where 70 is a timelike unit vector. The algebra of the vectors { 7 ^} is the same as 
that of Dirac’s 7 -matrices, but they now form an orthonormal basis for spacetime, 
rather than the four components of a single vector in an internal ‘spin-space’. 

By repeated multiplication we generate a basis of sixteen elements for the STA, 

1 (7m 1 {<7b ,i<Tk} (*7m) * 

1 scalar 4 vectors 6 bivectors 4 trivectors 1 pseudoscalar, 

where <r k = 7*70 and 


i = 70717273 = 0i<7 2 <7 3 . (2.2) 

The grade-4 element i is the highest-grade element in the space, and i 2 = — 1 . 
However, unlike the uninterpreted scalar imaginary used in conventional quantum 
mechanics, which we denote by j, it now has a definite geometric interpretation. The 
timelike bivectors { 07 } obey the same algebraic relations as the Pauli spin matrices, 
but in the STA they represent an orthonormal frame of ‘relative’ vectors in the 
Euclidean 3 -space relative to an observer whose 4-velocity is 70 . We shall take this 
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frame as our laboratory frame. To describe observations made in this inertial system 
we make use of the space-time split: if £ is a vector labelling an event in spacetime, 
this event is seen to occur at time t and position x where 

x^o = t + & = x '1o + (2.3) 

The ‘relative’ vector x is denoted in bold typeface to distinguish it from (proper) 
vectors. We denote the derivative with respect to spacetime position by V. This 
operator has the space-time split 70 V = dt + V, where V is the usual gradient 
operator in three-dimensional space, V = o^k- 

Two useful operations are reversion, denoted by a tilde, which reverses the order 
of 4 -vectors in an STA expression, and Hermitian conjugation, denoted by a dagger, 
which reverses the order of relative 3 -vectors in a non-relativistic expression. Angled 
brackets (tp)k are used to denote the grade -k part of the enclosed expression, with 
the subscript zero suppressed when discussing the scalar part. Our notation and 
conventions follow those in [ 8 , 10 ]. 

We shall need the STA representation of the Faraday tensor, in order to discuss 
the interaction of an anomalous magnetic moment with the electromagnetic field. 
The components may be assembled with the set of bivectors { 7 ^ 7 "} to form 
the Faraday bivector F, given by 

F = VAA = § 7 ^ 7 "^. (2.4) 

This has the space-time split F = E + iB, where E and B are the electric and 
magnetic fields in the 70 system. 

The even sub-algebra of the STA comprises a representation for the Dirac spinor 
■ 0 , allowing a reformulation of Dirac theory in the STA. This formalism and some 
applications to electron physics have been discussed at length in [ 11 ], which includes 
a preliminary summary of the results presented here. 

In this representation we may write [12] 

0 = (pe^R, (2.5) 

where p is a positive scale factor and R is a spacetime rotation (rotor) satisfying 
RR — 1. A direct method of translation between matrix-based approaches and the 
spacetime algebra of Hestenes is given in Ref. [10]. In the STA, the Dirac equation 
may be written as [7] 

V^zcr 3 — eAijj = m07o- (2-6) 

The components of the Dirac current assemble with the { 7 ^} vectors to give the 
conserved frame-free current J — ^> 70 $. The integral curves of J define a set of 
non-intersecting streamlines in spacetime, whose tangents are always future-pointing 
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timelike vectors. According to the standard Born interpretation of quantum mech¬ 
anics, J°(x,t ) gives the probability of locating the particle at position x and time 
t. The conservation of J implies that 

J°(xo,to)d 3 xo = (2.7) 

where (f 0 ,® 0 ) and (ii,®i) are joined by a streamline. The density J° thus flows 
along the flux tubes formed by adjacent streamlines, without ‘leaking’ out. For this 
reason, the streamlines are a useful tool for studying the flow of density in dynamic 
situations, independent of whether one accepts Bohm/de Broglie theory or any other 
particulate interpretation. It is in this spirit that we plot streamlines for the spin 
measurement in our simulations. The components (0|7/ys|0) of the spin current, 
like those of the Dirac current, assemble with the { 7 ^} vectors to give the frame-free 
spin current ps = ' 073 '*/’. 

The above discussion suggests that the role of the Dirac matrices in the theory is 
simply to define an orthonormal frame of vectors. It follows that claims of dynam¬ 
ical consequences based solely on their algebraic properties are questionable. One 
such example is the deduction that the measured velocity of an electron will always 
yield the speed of light (because the eigenvalues of the ‘velocity’ operators 707 * are 
±1). To the extent that we may sensibly define a particle velocity v in the Dirac 
theory, consistency with the laws of quantum mechanics requires that the velocity 
be in the direction of the Dirac current, i.e. J = pv. If not, particles would end 
up at positions with probabilities inconsistent with those derived from the Dirac 
(charge) current. Similar arguments may apply to spin measurements, where the 
anti-commutation properties of the Dirac matrices (or, non-relativistically, the Pauli 
matrices) are usually taken to imply that no two components of the spin may be 
known simultaneously. 

Our starting point is a form of the Dirac equation which describes an electrically- 
neutral spin-| particle of anomalous magnetic moment p. In the conventional rep¬ 
resentation this equation may be written as [13] 

Urd, - \^F. V ) | 0 (a 0 ) = m| 0 (®)). ( 2 . 8 ) 

The operators a 111 ' are defined in terms of the Dirac 7 -matrices by 

^ = ij[r,7l. (2.9) 

Translating equation ( 2 . 8 ) into the STA we find the equation 

V0icr 3 — — m07 O . (2.10) 

Note that there are some inconsistencies in the literature with the sign both of F 
and the term containing it in the modified Dirac equation. For example, compare 
Section 9.11 of Greiner [13], with equation (10.81) of Bjorken and Drell [14]. For 
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this reason, and to make contact with the non-relativistic treatment, we demonstrate 
that the choice made here is consistent with the Pauli equation in a non-relativistic 
reduction. Let us define two spinors ip± where 

ip±(x,t) = i(^±7o^7o)e i<r3m< - ( 2 - 11 ) 

The factor e ia3Tnt removes the contribution of the rest mass to the phase of ip± 
and, for s uffi ciently slow temporal variation in ip±, ensures that we are dealing with 
positive-energy solutions of (2.10). We identify the spinor ip 4 . with the Pauli spinor 
$, since it is an even element of the geometric algebra of Euclidean 3-space (the 
Pauli algebra [ 8 ]). Equation (2.10) then implies that 

2 mip- —V$i(j3 — fiEQicr 3, (2.12) 

where we have neglected the terms dtip-io'z and iiBip—cr 3 , which are small compared 
to 2 mip- in non-relativistic situations. Using this approximation for ip _, we find the 
equation of motion for the Pauli spinor: 

< 9^3 = —— (V 2 $ + nV(E^) - fiEV<S> - n 2 E 2 ^) - iiB<S>a 3 . 

2m V (2.13) 

The terms involving E are relativistic corrections to the Pauli equation arising from 
the motion of the magnetic moment through the electric field. In the limit c —► 00 
these terms vanish, leaving the STA form of the Pauli equation. In the STA, the 
parametrisation of $ used by Dewdney et al. (equation (1.1)) takes the transparent 
form 

$ = ^i/2 e *V 3 ^/2 e *Vi^/2 e i<7 3 V>/2 j (2.14) 

where the rotor term is precisely that needed to describe a rotation in terms of Euler 
angles. 

Using equation (2.12) it is readily shown [3] that a space-time split of the Dirac 
current gives a 3 -current J where 

mJ = — V ($ — |V-(3>i<7 3 $t), (2.15) 

and we have taken E = 0. The overstars in this equation denote the action of the 
gradient operator. The first term is just the STA equivalent of the current employed 
by Dewdney et al. The second term can be written as |Vx($$ 1 "s), which includes 
the curl of the spin 3-current <3>$ t s, obtained as an approximation to the Dirac spin 
current on splitting in the 70 system. This spin 3 -current has components given 
by equation (1.4). Hence the current initially employed by Bohm et al. [ 6 ] and 
more recently by Dewdney et al. is inconsistent with the non-relativistic reduction 
of the Dirac theory, although, as already noted, this does not affect their qualitative 
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conclusions. This fact seems to be widely known [4] but is still ignored in many 
causal treatments of the Pauli equation. The Dirac current arises naturally as the 
conserved current conjugate to 17(1) gauge transformations of the Dirac field, and 
hence appears as the source term in the QED Lagrangian. It is thus natural to adopt 
this definition of the current in any causal approach to Dirac theory. 


3 A model for the spin measurement 

Suppose we wish to measure the 2 component of the spin of an uncharged spin-| 
particle. Let the mass of this particle be m and the magnetic moment //,. We observe 
the deflection of the particle after passage through an inhomogeneous magnetic field 
directed along our 2 coordinate axis, and so infer the component of spin along the 
field direction. We will take our laboratory system to be the y 0 system. 

Let the initial state of the particle be a positive-energy plane-wave solution of 
the free-particle Dirac equation. We write this state in the form ij> = ip 0 e~ t<T3P ' x , 
where p is the momentum and the constant spinor r (po satisfies 

Pi >0 = rn*l> o7o- (3.1) 

An explicit form for ifio is (see e.g. Ref. [15] ) 

-ipo = e“ /2 $, (3.2) 

where $ is an arbitrary Pauli spinor and 

P7o = E + p = me u . (3.3) 

E and p are the energy and momentum of the particle in the q 0 system. We follow 
Dewdney et al. [ 1 ] by assuming that, at time t = 0, the particle receives an impulse 
from a magnetic field gradient, F = B’zia 3 8(t ), where B' is the field gradient along 
the ^-direction. Clearly this field is not a physically realisable magnetic field, since it 
has non-vanishing divergence. This could be rectified by the addition of components 
along the orthogonal directions with suitably chosen spatial gradients. However, 
we believe that such components would have little effect on the motion along the 2 
direction, and so for simplicity we ignore them. Around this impulse, the equation 
of motion for the system is given by the split of equation ( 2 . 10 ) in the 70 system, 
which approximates to 


d t ipia 3 = Apz8(t)’y 3 tp'y 3 , ( 3 . 4 ) 

with Ap = pB'. It follows that 

d t (i> ± 73^73) = T^pz8(t)(rp ± 73^73)^3- (3.5) 
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We may summarise the effect of the impulse by 

(if> ± 73^73)t=o + = (tp ± 73'073)t=o_e =FApz * (T3 , (3.6) 

where t = 0 _ is just before the impulse and t = 0+ just after. If we decompose the 
spinor i \)o into 

= KV’o - 73^073)5 ^ = |(^o + 73^073) ( 3 - 7 ) 

we see that, immediately after the shock, xp is given by 

</, = ^t e «r 8 (p-* + Apz) + 0l e tV 3 (p-* - Apz)_ (3. 8 ) 

The spatial dependence of ^ is now appropriate to two different values of the 3- 
momentum p 1 and fd, where 

fd = p + Apa^, fd = p — Apa^. (3.9) 


For t > 0 we have F = 0 , so that ip propagates as a linear combination of four 
plane-wave states 


<P = + rple icr3 P^ x + e^’*, 

(3.10) 

where 


p T 7o = E t + p T 

(3.11) 

"“ft, 

1 

II 

7a, 

(3.12) 

and 


E T = (m 2 -fp^ 2 ) 1 / 2 , -Ed = (m 2 + _p-*' 2 ) 1 / 2 . 

(3.13) 

with similar expressions for pb Each term in equation (3.10) must separately satisfy 

the free-particle Dirac equation, so that 


^3 + 

II 

+-* 

(3.14) 

— p^ipl. = m?/d 70 , 

(3.15) 

which are satisfied together with the condition 


<-i 

+ 

II 

7§- 

(3.16) 

These equations hold also for ip 1 , ip+ and ip}_. From these conditions we 

find that 

ip+ = 2£;T (p T 7o^ T + m7<W> T 7o), 

(3.17) 

ipi = 2£ , T (p T 7 o</> T m7oV> T 7o), 

(3.18) 
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with the expressions for \p± following similarly. 

We see that the elfect of the impulse on a monochromatic plane wave is to split 
it into four plane waves, each with distinct phase velocity. The first two terms in 
equation (3.10) have phase velocities in the directions ±(p + Apa 3 ) respectively, 
whilst the last two are in the directions ±(p — Apa 3 ). A particle incident on the 
apparatus could be modelled as a wavepacket, formed by superposing the incident 
plane-wave states above. This packet splits into four packets as a result of the B- 
field gradient. The two transmitted packets, which are composed of positive-energy 
plane-waves, move apart along the field direction until eventually they no longer 
overlap significantly, at which point detection of the particle at a screen would give 
an unambiguous result for the measurement. The two reflected wavepackets are com¬ 
posed of negative-energy plane-waves. The appearance of antiparticle states must 
ultimately be attributed to pair production in the magnetic field, and is significant 
only in large magnetic field gradients. A rigorous treatment of this phenomenon 
requires a multiparticle approach, which we shall not attempt here. 

We now verify that in the non-relativistic limit, this causal approach yields the 
same probabilities for a ‘spin-up’ and ‘spin-down’ result as conventional operator 
methods. In this limit, in which we neglect pair production, the negative energy 
terms in equation (3.10) are negligible and so we take ip\ = xft and ip\. = tpK If 
the initial state is a wavepacket composed of positive-energy plane-wave states, each 
having a weight g(p), then after passage through the magnetic field we have 

#M) = rf>u + t/> d = [ d 3 pg(p ) (Ve-^'* + , 

V J (3.19) 

where ip u denotes the packet moving up the field gradient and ipd the downwards 
moving packet. For sufficiently large t the overlap of these packets is small and may 
be neglected in calculating the density J°. In this approximation J° is given by 

J° = + ipdlo^pd- (3.20) 

The probability P u of obtaining the ‘spin-up’ result is simply equal to the probability 
of finding the particle in the region z > 0, which is given by integrating J° over this 
region. The contribution from ipd may be neglected since this packet is localised in 
the region z < 0. On performing the spatial integration we find that 

P u oc J d 3 p#(p) 2 (0 T 7o^ T )-7o, (3.21) 

with a similar expression for the probability Pd of the ‘spin-down’ outcome. Using 
equations (3.2) and (3.7), it is straightforward to show that 

(^ T 7 o^ T )- 7 o = 2 &$(E(p)/m + s-cr 3 ), (3.22) 

(V> i 7oV’ 1 )-7o = 2$ t $(£(p)/m - s-cr 3 ), (3.23) 
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where = $(j 3 $h If we take the Pauli spinor $ to be of the form 

$ = a — f3icr 2 , (3.24) 

and set Emm for all momentum components of the non-relativistic wavepackets, 
we recover the expected result that the ratio P u /Pd of the probabilities of the two 
outcomes is equal to |o| 2 /|/3| 2 . 

4 Wavepacket Simulations 

For computational simplicity we take the incident particle to be localised along the 
field direction only, and with no momentum components transverse to the field. This 
two-dimensional (one 3-space and the time coordinate) setup was used by Dewdney et 
al. [1] and is sufficient to demonstrate the essential features of the measurement 
process. The main deficiency of this model is that all four packets have group 
velocities along the field direction. 

We take ip 0 to be of the form (3.2), with u = ua 3 and $ given by equation (3.24). 
The parameters a and /? are real scalars, where a — 1,(3 = 0 corresponds to the 
‘spin-up 5 and a = 0, /? = 1 to the ‘spin-down 5 eigenstates of conventional operator- 
based formulations. These plane wave spinors are superposed numerically to form 
a Gaussian packet with standard deviation 1 x 10 _24 kgms _1 in momentum space, 
giving a half-width of 1 x 10 -lo m in real space. We take the particle’s mass to be 1 x 
10 -27 kg. After the impulse, the future evolution of ip is found from equation (3.10). 
From this information we can analyse the behaviour of the relative spin-vector 5A70, 
the probability density J Q = J- 70, and the streamlines for initial values of $. 

In Figure 1 we plot the evolution of Jo, of s A 70, and of streamlines starting 
from various positions in the initial packet whose spin-vector initially points along 
the ( 7 i direction ($ = 1 — ia 2 )• The impulse imparts A p = 1 x 10“ 23 kgms _1 . 
The packet splits symmetrically into two equal sized packets and the streamlines 
bifurcate at the origin. The plot of the spin vector shows that immediately after the 
impulse, the spins are disordered because of the infinite rate of spin precession in 
the impulsive field. However, they sort themselves out in around 10 -14 s such that, 
in the upward moving packet, the spin vector points along the field direction, whilst 
in the downward moving packet they are aligned in the opposite sense. Thus we 
recover the disjoint outcomes conventionally obtained by projection of the initial ip 
onto eigenstates of the spin operators. As Dewdney et al. found, the choice of which 
packet a streamline enters is determined by its starting position in the initial packet. 
If, for example, the streamline is initially above the bifurcation point, it enters the 
upward-moving packet and the spin evolves along this path such that, on arrival at 
a particle detector, the spin points along the field direction. We see that, for this 
choice of $, each outcome is equally likely, as is shown by the symmetrical splitting 
of the packet. 
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Figure 1: Splitting of a wavepacket with $ = 1 — icr 2 . Figure la show the probability 
density, J$, at t = 0,1.3,2.6,3.9 X 10 _14 ,s with t increasing up the figure, and Fig¬ 
ure lb the streamlines in the (t,z) plane. The evolution of the spin vector projected 
in the (x,z) plane is shown in Figure lc. 








z (Angstroms) 


z / A 

( 2 a) 


time / 10" u s 

( 2 b) 



( 2 c) 


Figure 2 : The asymmetric splitting of a wavepacket with $ = 1.618 — icr 2 . The 
bifurcation point is now below z = 0 so that more streamlines enter the upward 
moving packet. Figure 2c shows the evolution of s A 70 projected in the (x,z) plane. 






A calculation of the relative amplitudes of ip]. and ip. j_ shows that the negative 
energy packets are negligible for our chosen value of Ap, which is consistent with 
the observation that A p <C m implying negligible pair production. 

In Figure 2 we plot the evolution for unequal mixtures of ‘spin-up’ and ‘spin- 
down’ in the initial spinor. The impulse is as before. We now see an asymmetry in the 
splitting of the wavepacket, reflected in the lowering of the bifurcation point. More 
streamlines enter the upward moving packet and the density J 0 is correspondingly 
greater in this packet. We can calculate the probability for each of the two disjoint 
outcomes by numerically integrating J 0 over the spatial extent of each packet, giving 
a result in exact agreement with the results of the usual ‘projection onto eigenstates’ 
approach, as expected from our previous proof. 

As a final illustration of our approach, we consider a large field gradient such that 
A p is comparable to the mass of the particle. We suppose that the initial packet 
already has its spin aligned along the field gradient ($ = 1). For small field gradients 
we find that there is no splitting of the packet, all of the streamlines entering an 
upward-moving packet. However, for sufficiently large field gradients we find that a 
second packet is created which is composed of negative-energy states. In Figure 3 we 
plot the effects of an impulse Ap = 1 x 10 -18 kg m s -1 . Both packets have spin aligned 
with the field gradient, but the packet composed of negative-energy states is deflected 
down the field gradient. The negative-energy states therefore behave as if they 
have a ratio of magnetic moment to mans opposite in sign to their positive-energy 
counterparts. In a field-theoretic treatment the appearance of negative-energy states 
would be attributed to pair production, with the antiparticles having the opposite 
sign of magnetic moment but positive mass and giving the opposite sign for their 
ratio, as we have observed. 


5 Conclusions 

Our aim in this letter has been to present a fully relativistic, causal account of the 
operation of a Stern-Gerlach apparatus. This analysis extends the calculations of 
Dewdney et al. [1] to the relativistic domain, but with the important diff erence that 
our model does not use any specific physical interpretation of quantum mechanics. 
Our plots of the streamlines and spin vectors are in qualitative agreement with those 
of Dewdney et al. in non-relativistic situations, despite their use of a current that is 
inconsistent with the Dirac current in such cases. However, our interpretation of the 
streamlines is very different from that of Bohm/de Broglie theory, where streamlines 
are taken to be the spacetime trajectories of a statistical ensemble of particles. Our 
use of streamlines is solely as a visual aid to monitor the flow of the density J°. In 
the example considered here, the streamlines allow us to follow the flow of density 
from the initial wavepacket, through the region of interaction with the magnetic 
field, and into the final wavepackets that emerge from the apparatus. A similar use 
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Figure 3: Creation of negative-energy states by a large field gradient, giving an im¬ 
pulse of Ap = 1 X 10“ 18 A#tos“ 1 . The initial packet is entirely spin-up ($ = 1) and 
the resulting spin vector field (not shown) points uniformly along the cr 3 direction. 
Figure 3a is the probability density, Jo, plotted at t — 0,3, 6 ,9 x 10 _ 19 <s ; and Fig¬ 
ure 3b shows the streamlines in the (t, z) plane. The small packet moving along the 
—(73 direction is composed of negative-energy states. 


of streamlines in tunnelling simulations have allowed us to make novel predictions 
for the duration of electron and photon tunnelling processes. These results will be 
presented in a forthcoming paper. 

Our simulations make it clear that the Stern-Gerlach apparatus acts primarily 
as a polariser. It is the fact the dependence of the relative densities of the final state 
wavepackets on the incident spin that allows the experimenter to infer information 
about the incident spin. We could extend our simulations to include successive spin 
measurements, for example by directing particles in one output channel of an appar¬ 
atus with the field along <j 3 , into another apparatus with field along the < 7 i direction. 
This would not present any difficulty; we would simply find that the packet selec¬ 
ted from the first apparatus would divide symmetrically at the second apparatus. 
Having developed a relativistic model for a single particle spin measurement, we 
intend to extend the treatment to describe correlated measurements made on pairs 
of particles with spacelike separation. This should shed much light on the issues of 
non-locality that EPR and other experiments raise. 
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